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Abstract 

Using density-functional calculations we study the structure and polarization response of tetrag- 
onal PbTi03, BaTi03 and SrTi03 in a strain regime that is previously overlooked. Different 
from common expectations, we find that the polarizations in all three substances saturate at large 
strains, demonstrating a universal phenomenon. The saturation is shown to originate from an 
unusual and strong electron-ion correlation that leads to cancellation between electronic and ionic 
polarizations. Our results shed new insight on the polarization properties, and reveal the existence 
of a fundamental limit to the strain- induced polarization enhancement. 

PACS numbers: 77.22.Ej, 77.65-j, 77.80.-e 



Response of electrical polarization to external strains in infinite solids is a subject of inter- 
est from both fundamental and technological points of view.{l] Fundamentally, the process 
is governed by electrons, ions and their complex interactions, forming an important class of 
collective phenomena. Seeking direct physics and/or mechanism underlying the responses 
is of key relevance. Previously, Ederer and Spaldin showed that the strain dependence 
of polarization in different substances (BaTiC>3, BiFe03, LiNbOs, and PbTiOa) could be 
understood by the piezoelectric and elastic constants of the unstrained materials. [2| Techno- 
logically, inplane compressive strains were able to dramatically raise the critical transition 
temperature in BaTi0 3 (Reffl) and to turn incipient SrTiC>3 into strongly ferroelectric 
These important studies focus on small strains near equilibrium. 

Nowadays, advance in experimental fabrication of various high-quality ferroelectric(FE) 
nanostructures^] (such as thin films, nanowires, and dots) opens an opportunity for studying 
polarization response in a new regime with medium and/or large strains that have not been 
previously investigated {2I, [3, Q]. Unlike bulks, FE wires or dots allow exceedingly large 
external strains without formation of dislocation, since stress can be effectively relaxed due to 

Q 

finite lateral size. .6] For instance, large compressive strains can be realized in FE nanorods[7] 
by compressing them along the lateral directions. Even for thin films, strains of more than 
3% were shown possible in LaA103/SrTi03.[8( The new strain regime is interesting and may 
reveal many new physics in the sense that the charge density is heavily deformed by large 
strains, making electrons behave more collectively in some spatial region, but less in another. 
As a result, answers to questions (1) what determine FE off-center displacements in the new 
strain regime, and (2) how piezoelectric response, as a collective phenomenon, adjusts itself 
to large strains, remain unknown and are important. Furthermore, based on the widely 
accepted conclusion that polarization is to be enhanced as compressive strain increases 2, 

Q| , scientists have been long wondering what is the upper limit that polarization is able 
to reach, provided that exceeding strain is possible in laboratory. This is an intriguing 
question since it defines a profound limit on our widespread effort of seeking enhanced 
polarization by strain. 3, J] 

Here by means of first-principle calculations we investigate the structure and polarization 
responses of FE perovskites in the new strain regime. Three paradigm materials [namely, 
PbTiC>3 (PT), BaTiC>3 (BT), and SrTiC>3 (ST)] are studied, aimed to obtain conclusions 
that are generally applicable. We find that polarizations in considered FEs show a universal 
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behavior, by exhibiting three distinct stages in their strain-induced responses (each stage 
producing dissimilar piezoelectric coefficients). Two out of these three stages are new and 
previously unknown. We further reveal that when the compressive inplane strain is high, 
the polarization will no longer increase, defying the common belief that polarization must 
increase with increasing compressive strain. In fact, our calculations predict the existence of 
a critical inplane strain, above which the polarization starts to saturate. Our results suggest 
that there is a fundamental constraint on polarization enhancement by mechanical strain. 

Three considered materials represent FEs of different magnitudes of polarization: PbTiC>3 
is strongly ferroelectric, BaTiC>3 is ferroelectric but with a weaker polarization, and SrTiC>3 
is not FE under zero strain. We study the structural phase of tetragonal symmetry 
( |a.! | = |a 2 1 =a, |a 3 |=c), in which materials can have polarizations along the [001] direction. 
Biaxial in-plane strain (in Voigt notation) is defined as rji=r]2=(a — ao)/ao, where ao is the 
unstrained in-plane lattice constant. 

Polarization in infinite periodic solids can be expressed in lattice-vector coordinate as 
P = ^(xi&i + %2 a 2 + X3 a 3)> where fl is the volume of a unit cell and Xi describes the po- 
larization along the direction of lattice vector aj. Displacement of ions and deformation of 
electron density both contribute to P. Computing of ionic contribution (labeled as Pi on ) is 
straightforward using point charges, while electronic contribution (denoted as P e i) is deter- 



mined using the modern theory of polarization |10l. as implemented in our mixed-basis 
method. More specifically, P e i is computed as the geometrical Berry phase of valence elec- 
tron states as P e] = i-^ J dk 1 V k2 ^»(k 1 , k 2 )| kl =k 2 , where 0(k 1; k 2 ) = lndet| < u mill \u nk2 > \ 
is the phase of the determinant formed by occupied Bloch wave functions u nk . When inplane 
strains are imposed, two factors give rise to the change in polarization: (1) the variation of 
Xi, and (2) the cell-structure modifications of Q and a^. The first factor is critical and is 



manifested 
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13] in the proper piezoelectric coefficients as e^ 7 = ^ ^-a iai where is 



the strain in Cartesian notation. The second factor contributes to the improper piezoelectric 
coefficient which is trivially related to the proper one by the magnitude of polarization as 
e^ ro = e^™ 7 + 5 a /?P 7 — S/3yP a . One thus notes that Xi is a critical quantity that is tech- 
nologically relevant in terms of obtaining high piezoelectric efficiency This study will be 
concerned with the strain dependence of x%- 

Technically, we use density functional theory within the local density approximation 
(LDA)Q to determine optimized cell shape and atomic positions by minimizing the total 
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energy. Ca.cu.ations are performed using pseudopotentialmethod with mixed-basis set.H 
The Troullier-Martins type of pseudopotentials are used|16|]; atomic configurations for gen- 
erating pseudopotentials, pseudo/all-electron matching radii, and accuracy checking were 
described elsewhere. 1^] For each in-plane strain, the out-of-plane c lattice constant and 
atomic positions are fully relaxed. Our LDA-calculated inplane lattice constants of un- 
strained bulks are a=3.88A for PT, a=3.95A for BT, and a=3.86A for ST, which all agree 
well with other existing calculations. 

Structural response: Fig{T](a) depicts the c/a ratios of the three considered materials 
under compressive strains. Let us examine the BaTiC>3 curve first. For this curve, (1) our 
calculations confirm the (expected) linear relation between c/a and a when inplane strain 
is small; (2) However, as the inplane a length is decreased to 3.82A [point A in Figfjja)], 
BaTiC>3 apparently enters a different state, as witnessed by the fact that the c/a rises more 
sharply afterwards; (3) Interestingly, when a reaches point B, the c/a becomes linear again. 
These results combine to make the c/a in the region between A and B appear to be notably 
enhanced. This enhanced c/a regime in FEs has not been emphasized before (since previous 
studies deal with small strains) — and provides an explanation for the (mysterious) enlarged 



electromechanical response recently observed in PZT films under very large electric fields 
because these fields drive the films into the new strain range between A and B. 

Contrast of three materials in FigJIJa) tells us that the above conclusions are applicable 
also for PbTiC>3 and SrTiC>3. But the c/a value and the location of point A reveal interesting 
differences. First, we predict that for a fixed a constant, BaTi03 exhibits a larger c/a than 
PbTiC>3, unlike the zero-strain case where PT has a much higher c/a ratio of ~1.04 (in 
LDA). Meanwhile, the c/a in SrTiC>3 is considerably lower than in BT or PT, which could 
be attributed to the small size of Sr. When strained, the big Ba and Pb atoms repel other 
atoms more forcefully, leading to large c lengths. This size difference is also reflected in the 
location of the critical A point (at which the c/a changes behavior). The a value of point A 
is remarkably similar for PT and BT, but not for ST. 

FigJT](b) shows atomic off-center displacements (in unit of c) in BaTiC>3 under differ- 
ent strains, obtained by placing the Ba atom at the cell origin and then determining the 
displacements of Ti, 01 and 02 atoms with respect to their high-symmetry locations in 
tetragonal structure. Labels of individual atoms are given in the inset of FigJT^b). Atomic 
displacements in FigOJb) explain why point A plays an important role in determining the 
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structural response. When a > (lower strains), the O displacements are small. But, as 
a is below a^, the oxygen off-center shifts increase dramatically with a markedly different 
slope. Furthermore, the Ti displacement is positive (out of phase with the O displacement) 
when a > a^, and then becomes negative (in phase with the O displacement) for a < cla- 

Given the complexity in atomic displacements as well as in the c/a ratio, one would desire 
to find a mechanism (if any) that may determine how, and in what magnitudes, atoms are 
displaced. For this purpose, we examine the bond lengths in strained BaTiOs, shown in 
Figfjjc). Two observations are ready: (1) The Ti-Ol distance correlates well with the c 
length, (2) Despite that a is significantly compressed [by ~0.5A in FigOJa)], the Ti-02 and 
Ti-03 distances remarkably maintain to be nearly constant (so are the Ba-Ol and Ba-02 
bondlengths). Our calculations thus reveal the existence of a simple (and preferred) rule 
that can quantitatively explain the structural response over a wide range of inplane strains, 
that is, FEs respond to the external strains in a manner that maintains the lengths of short 
Ti-0 and Ba-0 bonds. This rule is found true also for PT and ST, and may be easily verified 
in x-ray experiments. 

Nevertheless, FE perovskites do not maintain the bond lengths in a trivial way. In 
Figfjjd) we describe the strain dependence of the optimal cell volume. Though cell volume 
is often assumed to be constant under biaxial strains, obviously this is not the case for FEs. 
Instead, Figfjjd) displays features that are very suggestive in terms of understanding the 
structure and polarization responses. Under small strains, the cell volume first decreases. 
This decrease is most pronounced in SrTiOs. Then approximately at a a, the volume of 
the material increases quickly, and later at as, it starts to decrease again. One thus sees 
that between a a and a^ the enhanced c/a in Figfjja) — and the enlarged O displacements 
in FigQJb) — occur by the expansion of cell volume. 

Polarization response: Figj2](a) depicts the Xs polarization and the proper e^ piezo- 
electric coefficient in PbTiOs under varied compressive strain. Caution is made to ensure 
that we follow continuously the same polarization as the X3 crosses different branches, since 
polarization in periodic solids is a multi-valued quantity. At zero strain, PbTiOs shows a 
theoretical Xs of 0.63 that corresponds to a polarization of 66.7/xC/cm 2 , consistent with the 
experimental result 



19| of ~70/iC/cm 2 . Upon strain, the most notable result in Figj2](a) is 
that PbTiOs displays three stages (labeled as I, II, and III) in its polarization response. 
At stage I (when a is above 3.84A), the X3 polarization increases linearly, from which 
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the piezoelectric e^ coefficient is deduced as 637 juC/cm 2 . At stage II, the X3 ~ a relation 
remains linear, but shows a different slope from that of stage I. We numerically find that 
the e^[° for stage II is 898 /zC/cm 2 . Interestingly, when the inplane strain continuously 
transitions from stage I to II, we notice an abrupt rise in the polarization. This leads to the 
occurrence of a huge piezoelectric e^ coefficient of 3900 /xC/cm 2 near a = 3. 84A. 

When the inplane a constant is further decreased below 3.70A (which corresponds to a 
moderate 4.5% strain), another unusual phenomenon occurs in Figj2](a). That is, the in- 
crease in polarization dramatically slows down, and the Xs finally saturates at a value of 
1.05. Further increase of inplane strain no longer affects the X3 polarization. This satura- 
tion is surprising, since it is generally accepted that FE polarization always increases with 
increasing compressive strains. The saturated polarization also demonstrates the existence 
of a maximal limit up to which polarization can be enhanced by means of epitaxial strains. 

We have also performed polarization calculations for BaTiC>3 and SrTiOs; results are 
shown in FigfS^b). For incipient SrTiC>3, we find a threshold strain at r]!=l%, below 
which the pol arization remains null. A similar magnitude of threshold strain was reported 
previously. 2o| Results for BaTiC>3 and SrTiC>3 in Figj2]^b) reveal that (1) the existence of 
a three-stage polarization response and (2) the saturation of polarization apply to all three 
substances, hence demonstrating a general phenomenon. Also interestingly, we find that X3 
saturates at similar values in BaTiOs and in SrTiC>3, but differently in PbTiC>3. 

The saturation of X3 can n °t be naively related to the saturations of effective charges 
and/or atomic displacements. First, we see in Figfjjb) that the O displacements do not 
saturate, in fact. To illustrate if effective charge saturates, we compute the effective Z% 3 
charge of each atom in PbTiC>3 under strains, by finite difference Z* = — ^j. Under zero 
strain, the effective charges (denoted as Zq) are 3.65 (Pb), 5.76 (Ti), -4.90 (01), and - 
2.28 (02). With strain, we are interested in the change of the effective charges, namely 
| Z* | — | Zq ], given in Figj2](c). Note that although the polarization increases once we 
turn on compressive strain, the effective charge nevertheless decreases for all atoms. In 
other words, the | Z* | has its maximum at the equilibrium structure. It is also interesting 
to recognize that the | Z* | — | Zq | in Figj2](c) separate into two groups: one group 
consists of Ti and 01 atoms, while the other is formed by Pb and 02 atoms. The strain 
dependence of the effective charge is very similar within the same group, implying the 
occurrence of a chemical correlation. Meanwhile, we note that the Z* do not saturate — and 



6 



for instance, the Z* of Pb (also of 02) slightly increases when the a length is below 3.7A. 
Since 6x3 = ^2i{^Z*Ar ijZ + Z*S(Ar ijZ )} (where i is the index of atom), the saturation in 
X3 thus results from a balanced cancellation between strain-induced variations of Z* and 
atomic displacements. 

We now provide microscopic insight into why polarization undergoes a (puzzling) change 
from stage I to stage II, and why X3 saturates at high strains, by examining electron charge 
distributions (Figj3]). The electron distribution in equilibrium PbTiC>3 is as expected: charge 
transfer leads to electron deficiency near atomic cores (blue regions), while more charge 
accumulates in the interstitial area between atoms (red regions). As PT is strained to 
a=3.80A, an important change occurs, namely the overlapping charge largely disappears 
between Ti and 01 atoms. This reveals that the (considerably) weakened Ti-01 bond is the 
reason responsible for the transition from stage I to stage II. It also explains why the c/a is 
enhanced in the second stage. 

On the other hand, when PbTiOa is further strained to a=3.72A, electron charge near 
the 03 atomic core is heavily deformed, as a result of two factors that (1) more electrons are 
transferred into the interstitial between Ti and 03, and (2) the strong Coulomb repulsion 
from these interstitial charges distort the electrons near the 03 nuclei. In this circumstance, 
electrons between Ti and 03 behave more like strongly correlated particles. Further, since 
the electrons near the 03 nuclei couple directly with this ion, ionic displacement and electron 
deformation are thus expected to be correlated as well. This expectation is indeed confirmed 
by the electronic (P e i) and ionic (Pi on ) contributions to the total (P to t) polarization [FigfSJ^d)]. 
When a decreases from ao to 3. 72 A, P e i rises faster than Pj on declines, giving rise to an overall 
enhancement of the total polarization. However, when a is below 3.72A, the P e i and Pi on 
contributions show equal slopes but with opposite signs (due to the above- described electron- 
ion correlation), thus leading to the polarization saturation in stage III. Another interesting 
indication of the strong electron-ion correlation is manifested by the fine structure on how 
the Ti-03 bondlength depends on strain [see the inset of Figj2](d)]. We find that when X3 
starts to saturate near a=3.72A, the Ti-03 bondlength surprisingly begins to increase. This 
demonstrates the existence of a minimum length (B min ) by which the Ti-03 short bond can 
be compressed. Our predicted B min is 1.76A (PT), 1.78A (BT), and 1.79A (ST). Reaching 
this minimum length implies the onset of polarization saturation. 

In summary, density functional calculations on strained FEs reveal two unusual properties 
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of polarization: (1) the occurrence of stage II with distinct behaviors, and (2) the saturation 
of polarization upon higher strains. We further predict that the continuous transition from 
stage I to II may produce marked piezoelectric responses. Microscopic insights responsible 
for these behaviors are discussed. The enhanced piezo-response and c/a ratio in stage II is 
associated with the weakening of the Ti-Ol bonds. On the other hand, the electron-electron 
and electron/ion correlations lead to the polarization saturation. Our calculations reveal 
the existence of a profound limit up to which polarization can be enhanced by strain. The 
predicted saturation value of polarizations are larger than 125 /zC/cm 2 for three substances, 
suggesting that currently there is still room for polarization enhancement. This work was 
supported by the Office of Naval Research. We thank H. Krakauer and D. Vanderbilt for 
very useful discussions. 

Note added: After having completed the manuscript, we were brought to a recent 
study 2l| where the polarizations in highly strained Pb(Zr .2Ti .8)O3 films (with c/a «1.09) 
is surprisingly similar to that in relaxed PZT films (with c/a ^1.05). This, in spirit, adds 
further evidence supporting the validity of our results on polarization saturation. Never- 
theless, our key conclusions that (1) polarization saturates in all materials, thus being a 
universal principle; (2) the saturation origins from strong electron-ion correlation; (3) there 
is an unusual volume expansion and an enhanced c/a structural response between points A 
and B, have not been addressed previously}^ . 
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FIG. 1: Structural properties as a function of inplane lattice constant: (a) the c/a ratios (where 
points A and B are marked only for BT and ST, for the clarity of display), (b) atomic displacements 
in BaTiOa, (c) bond lengths in BaTiO"3 (where the c-axis length is also plotted for comparison), 
(d) the change of cell volume AV/Vo (where Vq is the equilibrium volume). Labels of atoms in 
tetragonal perovskite are given in the inset of (b). 
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FIG. 2: Dependencies of the following properties as a function of the inplane lattice constant: 

(a) the X3 component of polarization (using the left vertical axis) and piezoelectric coefficient 
(using the right vertical axis) in PDT1O3, (b) the X3 polarizations in BaTiOs and SrTiO-3, compared 
with that in PbTiO"3, (c) effective charges | Z* \ — \ Zq | in PDT1O3, and (d) electronic (P e i) and 
ionic (Pion) contributions to the total (Ptot) polarization in PbTiC>3. For each calculated point in 
(d), the A-site atom is always placed at the cell origin so that comparison of P e i at different lattice 
constants are meaningful. Inset of (d) shows the fine structure of Ti-02 and Ti-03 bond lengths 

(b) versus the inplane (a) constant. 
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a=3.88 A a=3.80 A a=3.72 A 

FIG. 3: (Color online) Charge-density difference Ap = p — po, where po is the sum of electron 
density of free atoms, in PbTi03 at three different inplane lattice constants of a=3.88, 3.80, and 
3.72A. The plotted density is on the plane that cuts through the Ti, 01, 02, and 03 atoms. 
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